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We study the effect of resonant tidal excitation of inertial modes in neutron stars during binary 
inspiral. For spin frequencies less than 100 Hz, the phase shift in the gravitational waveform associ- 
ated with the resonance is small and does not affect the matched filtering scheme for gravitational 
wave detection. For higher spin frequencies, the phase shift can become significant. Most of the 
resonances take place at orbital frequencies comparable to the spin frequency, and thus significant 
phase shift may occur only in the high-frequency band (hundreds of Hertz) of gravitational wave. 
The exception is a single odd-paity m = 1 mode, which can be resonantly excited for misaligned 
spin-orbit inclinations, and may occur in the low-frequency band (tens of Hertz) of gravitational 
wave and induce significant (3> 1 radian) phase shift. 
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I. INTRODUCTION 

Coalescing neutron star-neutron star (NS-NS) and 
neutron star-black hole (NS-BH) binaries are the most 
promising sources of gravitational waves (GWs) for 
ground-based detectors such as LIGO and VIRGO Jj. 
The last few minutes of the binary inspiral produce 
GWs with frequencies sweeping upward through the 10- 
1000 Hz range (LIGO's sensitivity band)01. Due to 
the expected low signal-to-noise ratios, accurate gravita- 
tional waveforms are required to serve as theoretical tem- 
plates that can be used to detect the GW signal from the 
noise and to extract binary parameters from the wave- 
form. 

In the early stage of the inspiral, with the GW fre- 
quencies between 10 Hz to a few hundred Hz, it is usu- 
ally thought that the NS can be treated as a point mass, 
and tidal effects are completely negligible. This is indeed 
the case for the "quasi-equilibrium" tides as the tidal in- 
teraction potential scales as a~^ (where a is the orbital 
separation) [^ . The situation is more complicated for the 
resonant tides: As two compact objects inspiral, the or- 
bit can momentarily come into resonance with the normal 
oscillation modes of the NS. By drawing energy from the 
orbital motion and resonantly exciting the modes, the 
rate of inspiral is modified, giving rise to a phase shift 
in the gravitational waveform. This problem was first 
studied [J, Q in the case of non-rotating NSs where the 
only modes that can be resonantly excited are g-modes 
(with typical mode frequencies ^ 100 Hz). It was found 
that the effect is small for typical NS parameters (mass 
M ~ IAMq and radius i? ~ 10 km) because the coupling 
between the g-mode and the tidal potential is weak. Ho & 
Lai (1999)|6j studied the effect of NS rotation, and found 
that the g-mode resonance can be strongly enhanced even 
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by a modest rotation (e.g., the phase shift in the wave- 
form A<I> reaches up to 0.1 radian for a spin frequency 
Vs ^ 100 Hz). They also found that for a rapidly rotat- 
ing NS [vs ^ 500 Hz), f-mode resonance becomes pos- 
sible (since the inertial-frame f-mode frequency can be 
significantly reduced by rotation) and produces a large 
phase shift. In addition, NS rotation gives rise to r-mode 
resonance whose effect is appreciable only for very rapid 
(near breakup) rotations. Recently, Flanagan & Racine 
(2006) l7j studied the gravitomagnetic resonant excitation 
of r-modes and and found that the post-Newtonian effect 
is more important than the Newtonian tidal effect (and 
the phase shift reaches ~ 0.1 radian for Vs ~ 100 Hz). 
Taken together, these previous studies suggest that for 
astrophysically most likely NS parameters (M ~ IAMq, 
i? ~ 10 km, Vs ^ 100 Hz), tidal resonances have a small 
effect on the gravitational waveform during binary inspi- 
ral. 

A rotating NS also supports a large number of Coriolis- 
force driven modes named inertial modes (also called ro- 
tational hybrid modes or generalized r-modes; see, e.g., 
Refs. [alajll3jll3)i of which r-mode is a member. The 
modes have frequencies of order the spin frequency, and 
as we show below, they couple more strongly to the (New- 
tonian) tidal potential than the r-modc. It is therefore 
important to investigate the effect of inertial mode reso- 
nances on the GW phase evolution during binary inspiral 
— this is the goal of this paper. 



II. BASIC EQUATIONS FOR TIDAL 
RESONANCE DURING BINARY INSPIRAL 

We consider a NS of mass M, radius R and spin fig 
in orbit with a companion M' (another NS or a black 
hole). The orbital radius a decreases in time due to GW 
emission, Jlorb is the orbital angular frequency. We allow 
for a general spin-orbit inclination angle Q (the angle 
between ils and the orbital angular momentum L). In 
the spherical coordinate system centered on M with the 



Z-axis along L, the gravitational potential produced by 
M' can be expanded in terms of spherical harmonics: 
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[Here the symbol (— )^ is zero if p is not an integer.] Only 
the I > 2 terms are relevant, and the dominant (I = 2) 
tidal potential has W2±2 = (Stt/IO)^/^. Since it is most 
convenient to describe oscillation modes relative to the 
spin axis, we need to express the tidal potential in terms 
of Yim{d,(j)), the spherical harmonic function defined in 
the corotating frame of the NS with the z-axis along Jl^. 
This is achieved by the relation 
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where the !?„/„/ is the Wigner P-function (e.g., 
Wybourne 1974), and (j)g ^ (j) + fl^t. 

The linear perturbation of the tidal potential on the 
NS is specified by the Lagrangian displacement, ^(r, i), 
of a fluid element from its unperturbed position. In the 
rotating frame of the NS, the equation of motion takes 
the form 
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^ + 20.x^ + C.€ = -Vt/, 
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where C is a self- adjoint operator (a function of the pres- 
sure and gravity perturbations) acting on ^ (see, e.g., 
Ref.[lJ|)- A free mode of frequency cua with ^£j(r,i) = 
|„(r) e-*'^°* ex e™'^-*'^°* satisfies 



cjli^ - 2iu;Ms X ^„ + C • ^„ = 0, 



(5) 



where {a} denotes the mode index, which includes the 
azimuthal "quantum" number m. We carry out phase 
space mode expansionjl^, uM 
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Using the orthogonality relation[i3| (^„,2ir2s x ^^,) 

(Wa + Wa')(^a'^a') = ' 
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and we have used the normalization (Cq,^^) = 1. Recall 
that in Eqs. (|7I9|I . the index a includes m. 

Now consider the excitation of a specific mode with 
inertial- frame frequency CTq = Wq -t- mfls by the potential 
component oc e^"" * — we call this {a, m')-resonance. In 
the following we shall adopt the convention m > and 
to' > 0. We first consider the case of CTq > (i.e., the 
mode is prograde with respect to the spin in the inertial 
frame, although in the rotating frame the mode can be 
either prograde or retrograde, corresponding to cJq > 
and uja < respectively). The resonance occurs when 



aa — m r^orb- 



(11) 



Note that the mode azimuthal index m = m' for aligned 
spin-orbit (O — 0), but in general the TO.'-th potential 
can excite a mode with a different m. Integrating Eq. Q 
across the resonance, we find that the post-resonance am- 
plitude is given by 



= dtfa^r. 
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where a, J7orb should be evaluated at the resonance radius 
Ga, at which a^, = m'[G{M + M')/al]^^'^ . The mode en- 
ergy in the inertial frame is 2ctq£q|cqP, where the factor 
2 accounts for the fact that for each to, > 0, CTq > reso- 
nant mode (recall that our convention is to. > 0), there is 
also an identical to < 0, cTq, < resonant mode, and they 
should be counted as the same mode. Thus the energy 
transfer to the mode during the resonance is given by 
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where we have defined Qa.im in units such that M — R — 
1. Except for differences in notation and sign convention, 
Eq. l(T!^ agrees with the expression derived in Ref. 3, 
where mode decomposition was not carried out rigorously 

jma. 

The phase shift in the gravitational waveform due to 
the resonant energy transfer is twice the orbital phase 
shift and is given by 
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where i?orb 



-GMM' /{2a) is the orbital energy, and 
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is the orbital decay timescale {q = M' /M is the mass 
ratio; all quantities should be evaluated at the resonance 
radius Oq). Keeping the leading /-term in Eq. (jl3|l . we 
find 
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where a^ = cJaiR^/GMf/^ and £„ = edR^/GMy/^. 

For modes with ctq < (i.e., retrograde with respect 
to spin in the inertial frame), the GW phase shift at the 
resonance — CTq = m'Qorh can be similarly derived and is 
given by 
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III. 



INERTIAL MODES AND TIDAL 
COUPLING COEFFICIENTS 



To fix notations, we first summarize the basic prop- 
erty of inertial modes for incompressible stars[llj. To 
order 0{i^s), the mode displacement vector ^q, e~'"* ex 



satisfies the equation 



V^, 



(18) 



where e^ is the unit vector along the z-axis (Jlj), q = 
2Q.s/(jJ, and w^^ = 6P/p is the (Eulerian) enthalpy per- 
turbation. For incompressible fluid, V • ^^^ = 0, we have 



vV- 
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= 0. 
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This equation can be solved in an ellipsoidal coordinates 
(xi,a;2), giving 2p oc Pp{xi)Pp{x2)- For a given set of 
{m,j), there are (j — m) eigenvalues uja/^s- The modes 
with even (j — m) have even parity with respect to the 
equator, while those with odd {j — m) have odd parity. In 
this notation, "pure" r-modes correspond to those with 
j—m = 1 and have frequencies uJa — —2Q.s/{\+m). Such 
characterization of inertial modes can be generalized to 
stellar models which are compressible. 



A. Even-Parity Modes 



The m = 2, even-parity [j 
can be excited by the I = to' 



4, 6, • • •) inertial modes 
: 2 tidal potential. The 



frequencies of some modes are listed in Table I. Consider 
the two TO = 2, j — 4: inertial modes. For uniform stellar 
models, both have eigenfunctions 
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where fi = q ^ = a;/(2f2s), and {'nj,4>,z) are cylindrical 
coordinates. A direct calculation shows that 



Qc.,22 = j<fxpC-V{r^Y22) 

dip 



= /d^xp^— liu— - + 2^^ =0. (21) 
1 + q \ w 



Thus, to order 0{Q.s): the tidal coupling coefficient van- 
ishes. The first non-zero contribution arises in the order 
0{Vil). Integration by part of equation yields. 



Qa,22 = ld\5plr''Y22(e,4>) + ir^Y22pC-dS, 



(22) 



where Spa is the Eulerian density perturbation and 6pa — 
[ixf^p^ /T iP)il} with Fi being the adiabatic index. For- 
mally, Fi = oo in an incompressible model. The first 
term on the right-hand-side of Eq. H22I) scales as fig . The 
Lagrangian pressure perturbation vanishes at the surface, 
or 5P + ^ ■ VP = 0. So the surface radial displace- 
ment £^a,r satisfies ^q = {uj'^/g)ip{r = i?) oc O^ (where 
g = GM/R^). As a result, the surface integral on the 
right-hand-side also comes in at the order of 0{fl^). 

To accurately calculate tidal coupling to order 0{fl^), 
it will be necessary to compute the mode eigenfunction 
to order 0{fl^), including rotational distortion to the hy- 
drostatic structure and the centrifugal force in the per- 
turbation equation. This is beyond the scope of our pa- 
per. To estimate the tidal coupling coefficient, we con- 
sider a compressible, but uniform stellar model. The 
Eulerian density perturbation dpa is related to -0 by 
Spa = {ui'^p'^/TiP)ijj, and we fix Fi to a finite value. 
This 6 Pa is then used in Eq. it^ to obtain Qa,22- The 
results are given in Table I. 

For NS models with nonuniform density profile, no an- 
alytical solution for the inertial mode eigenfunction is 
generally possible. Wu 11] showed that if the star has a 
power-law density profile, p oc (i?^ — r'^)0 ^ equation l(T^ 
is separable in the ellipsoidal coordinates and one can 
easily obtain eigenfunctions that are accurate to O(fls)- 
However, the Eulerian density perturbation caused by 
inertial modes remains small [Spa/p — {uj'^/c^)ip ~ 
ujl/{GM/R^)V^tp < V • |„], we expect that Qa,22 is of 
order fl^, similar to that for the uniform density model 
(/3 — 0). Table I gives the mode eigenfrequencies and 
tidal coupfing coefficients calculated using the eigenfunc- 
tions from 11] and setting Fi = dlnP/dlnp. Correc- 
tions of order unity may arise when one is able to obtain 
eigenfunctions accurate to 0{ill). 



TABLE I: Frequencies and Tidal Coupling Coefficients for TABLE IL Same as Table |I] but for m — 1 inertial modes 
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0" 




:0.5 
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1.0 
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-1.2319 
0.2319 


0.084 
0.031 


-1.1607 
0.3830 


0.0027 
0.0031 


-1.1224 
0.4860 


0.013 
0.019 



6 -1.6434 0.018 -1.5820 0.0023 -1.5415 0.0022 

-0.8842 0.018 -0.8726 0.0077 -0.8671 0.0059 

0.1018 0.006 0.1780 0.0036 0.2364 0.0031 

1.0926 0.017 1.1814 0.021 1.2408 0.018 

"The stellar density profile is p oc {B? — r^)" . 

^u} is the mode frequency in the rotating frame, as in ^ ex; 

"Qc<,22 = ^d^xSpl (r^Yaz) = Qtli, with C, = C4i?VGM)i/2. 
Qa,22 is calculated with normalization id^xp^ • $ = 1 and M = 
R = 1. For the uniform-density model, Q is calculated fixing Fi = 1 
(otherwise the result is zero), while for other models we use Fi = 
din P/ din p. 



B. Odd-Parity Modes 

For an inclined orbit (0 ^ 0), the rn = 1, odd-parity 
(j = 2, 4, • • •) modes can also be excited by the I — m! ^ 
2 tidal potential (see §IV). The relevant tidal coupling 
coefficient is 



Q — Qa.,21 



d^xSplr^Y2i, 
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with F21 = "(15/87r)i/2sin6'cos6'e 

The m = 1, j — 2 mode is a r-mode. For a uniform 
incompressible star, the mode frequency is given by[l 
m, to order 0{nl), 



ujc/^s = -1 - 3^2^/4, Ua/VLs 
The tidal couphng strength isQ 
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For higher order (m =1, j = 4, 6, • ■ •) modes, eigen- 
functions including f2^ corrections are not available, 
and our results for Qa.21 are estimates using 5pa — 
{bj'^p^ /TiP)^, with a -0 that is accurate to 0{Vts). The 
mode frequencies and coupling coefficients are listed in 
Table II for the j = 4, 6 modes. 



IV. EFFECT OF INERTIAL MODE 
RESONANCE 

The result of §111 (see Table I) shows that the cou- 
pling coefficient of an inertial mode to the I — m — 2 
tidal potential has the form Qa,22 = Q^l, with Q ^ 0.1, 
and Cls = flsiR^/GMy/^. For most modes, Q is hkely 





13- 


= 


/3 = 


0.5 


13 = 


1.0 


j 


Lu/ils 


Q" 


Lo/^s 


Q 


Lo/Qs 


Q 


4 


-1.7080 


0.023 
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0.005 
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-0.6511 


0.009 
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0.012 
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0.052 



6 -1.8617 
-1.3061 
-0.4404 
0.5373 
1.4042 



0.005 
0.005 
0.030 
0.075 
0.008 



-1.8308 0.001 

-1.2941 0.006 

-0.4833 0.017 

0.6175 0.022 

1.4431 0.016 



-1.8088 0.001 

-1.2866 0.005 

-0.4126 0.001 

0.6726 0.019 

1.4688 0.011 



,„2i=Jd^xSp*^{r2Y2i) = Qnl 



to be significantly smaller. Given that our result for Q 
should only be considered as an order-of-magnitude es- 
timate (since we did not include fi^ correction to the 
mode eigenfunction), in the following we will scale our 
equations using Q ~ 0.1. Similar consideration applies 
to the tidal coupling of odd-parity modes (Table II). 

An m = 2 inertial mode is resonantly excited by the 
m' = 2 tide when aa = 0Ja + 2^}^ = 2J7orb- We only need 
to keep the Z = 2 tidal potential, with 1^22 = (Stt/IO)^/^, 
Vil^ I = cos4(e/2). Thus the GW phase shift due to the 
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where M1.4 = A//(1.4Mo), Rw ^ R/{10 km). Note that 
n,/ie^ujo^) ^ 1, and Cls = (i^./2170 Uz){RlJMi,iy/^ 
(where i^s is the spin frequency). Thus, for Vs ^ 300 Hz, 
^'I'gw ^ 0.05 radian (for i? ^ 10 km) — such a phase 
shift is too small to affect GW detection. 

Compared to the to = 2, j = 3 r-mode (often called 
I = m = 2 r-mode) resonance studied in Ref. Q, the 
phase shift associated with the inertial mode resonance 
(Eq. n^ is larger by a factor of (a/R)'^ oc Qs , because 
the r-mode can only be excited by the 171 — 2,1 — 3 
(Newtonian) tidal potential. 

For an inclined orbit (8 ^ 0), the I — m' = 2 tidal 
potential oc Y22{0L,4'L)e~^^'^ has a component propor- 
tional to I21 {&, 0), and thus can excite m = 1, odd-parity 
modes at the resonance \aa\ = 2r2oib- For ctq > (i.e., 
the mode is prograde with respect to the spin in the in- 
ertial frame), from Eq. H16|) with m = 1 and m' = 2, and 

I?^2^ = 2cos3(e/2) sin(e/2), we find that the GW phase 
shift associated with the resonance is 
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For modes with \aa\ ~ ^s (such as the m — 1, j — A 
modes), Eq. I|27|l has the same behavior as Eq. (|26|1 for 
the excitation of even-parity modes. The exception is 
the 171 = 1, j = 2 r-mode, which has (Ja/fts = — (3/4)rjg. 
Since this mode is retrograde with respect to spin in the 
inertial frame, we use Eq. H17(l to find 
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Thus the phase shift is always very significant for © not 
too close to 0. However, the resonant condition 2rJorb = 
|(Tq.| implies that the GW frequency at the resonance is 



.Gw = J^^-. = 10.2^' 



12. ( "'^ V 
1.4 V400 Hz/ 



Hz. 
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ri/3, 



Thus, only when i^s > 400Afj^(j /Rio Hz can the reso- 
nance occur in the LIGO band. 



V. DISCUSSION 

The result in this paper indicates that for neutron 
stars with spin frequencies Vs ^ 100 Hz, resonant exci- 
tations of inertial modes during binary inspiral produce 
negligible phase shift A$gw in the gravitational wave- 
form. At such spin frequencies, only the g-mode reso- 
nance 61 and the gravitomagnetic resonance of the m = 2 
r-mode |3| may produce A$gw of order 0.1 radian — 
such a phase shift probably does not affect the matched 
filtering method for detecting gravitational waves. Of 
course, one should keep in mind that the above num- 
ber is for "canonical" neutron stars with M = 1.4Mq, 
R = 10 km. The phase shift depends strongly on the 
neutron star parameters: For the g-mode resonance, 
A$GW Of R^-^M'^-^Vq^ (where j^gw is the gravita- 
tional wave frequency at resonance, and it depends non- 
linearly on the spin frequency); for the gravitomagnetic 



r-mode resonance, A<i>GW oc J?^M~-'^''/^i^Q^ (where i^gw 
of the same order as the spin frequency) . Thus the phase 
shift is larger for neutron stars with larger radii. 

For higher spin frequencies {vs ^ a few x 100 Hz), 
various resonances becomes important, including the f- 
mode resonances studied in Q and the inertial mode 
resonances studied in this paper. While most of these 
resonances occur at the high-frequency band of the grav- 
itational waves (I'GW is of order the spin frequency, and 
thus > a few x 100 Hz), it is of interest to note that the 
excitation of the m = 1 r-mode (the m = 1, j = 2 mode 
in our convention; see §HLB) occurs at the low-frequency 
band (tens of Hertz) of LIGO (see Eq. j^Hl) and produces 
very large phase shift (see Eq. [2H1)- While neutron stars 
with h's ^ a few x 100 Hz are found within binaries with 
a low-mass companion (e.g., white dwarfs), it is not clear 
that they are produced in compact binaries with another 
neutron star or black hole companion. If such systems 
do exist, then the detection the m = 1 r-mode resonance 
should be straightforward, and would lead to a clean con- 
straint on neutron star parameters (radius and spin). 

Finally, it is worth mentioning some caveats of the 
present study. Our calculations of the tidal coupling co- 
efficients for most inertial modes are approximate, since 
we did not take into account of the 0{Qf) correction 
to the mode eigenfunction. More importantly, we have 
neglected buoyancy in the stellar models. The Brunt- 
Vasala frequencyinside a cold neutron star could be as 
large as 100 Hz [a,|23- While the finite buoyancy has a 
negligible effect on the r-mode, it could affect the other 
inertial modes in an appreciable way. It is worthwhile 
to study the property (particularly the tidal coupling 
strength) of the inertial modes in nonisentropic neutron 
star models (see |21| and references therein). 
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